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ABSTRACT 

Episodic jets are usually observed in the intermediate state of black hole transients during 
their X-ray outbursts. Here we report the discovery of a strong positive correlation between 
the peak radio power of the episodic jet Pj et and the corresponding peak X-ray luminosity L x 
of the soft state (in Eddington units) in a complete sample of the outbursts of black hole tran¬ 
sients observed during the RXTE era of which data are available, which follows the relation 
log Pj e t = (2.2 ± 0.3) + (1.6 ± 0.2) x logL x . The transient ultra-luminous X-ray source in 
M31 and HLX- 1 in EXO 243-49 fall on the relation if they contain stellar mass black hole and 
either stellar mass black hole or intermediate mass black hole, respectively. Besides, a signif¬ 
icant correlation between the peak power of the episodic jet and the rate-of-increase of the 
X-ray luminosity dL x /dt during the rising phase of those outbursts is also found, following 
log P jet = (2.0 ±0.4) + (0.7 ±0.2) x logdi x /df. In GX 339-4 and H 1743-322 in which 
data for two outbursts are available, measurements of the peak radio power of the episodic 
jet and the X-ray peak luminosity (and its rate-of-change) shows similar positive correlations 
between outbursts, which demonstrate the dominant role of accretion over black hole spin 
in generating episodic jet power. On the other hand, no significant difference is seen among 
the systems with different measured black hole spin in current sample. This implies that the 
power of the episodic jet is strongly affected by non-stationary accretion instead of black hole 
spin characterized primarily by the rate-of-change of the mass accretion rate. 
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1 INTRODUCTION 

Radio jet has been observed in Galactic black hole X-ray binaries 
and neutron star X-ray binaries, in transient ultra-luminous X-ray 
binaries in nearby galaxies which are suspected to contain stellar 
mass compact objects or intermediate mass black holes (IMBHs), 
and in active galactic nucleus (AGN) which contains supermassive 
black holes (SMBHs). Accretion time scale in black hole X-ray bi¬ 
naries is much shorter than that of the AGNs, which allows us to 
study black hole accretion on longer accretion time scales unacces- 
sible in AGNs using current observations. For instance, we are able 
to observe distinct X-ray spectral states and spectral transitions in 
single black hole X-ray binaries with X-ray all sky monitoring or 
pointed observati ons. Dis tinct spectral states are probably exist in 
AGNs l lKording et alj2006h . but complete spectral state transitions 
have not been seen due to limited observing time scale compared 
with accretion time scales of SMBHs. Black hole X-ray binaries are 
therefore unique sources to probe accretion physics such as spectral 
state transitions. 

Most black hole X-ray binaries (BHXBs) are transients 
jLi li et al.ll2007 j). During outbursts, they will usually go through 
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several distinct X-ray spectral states, such as the hard state, in- 
termedi ate state (or ste ep power-la w state ) and the^high soft state 
(see iRemilliad & McClintockll2006l : iDone et ai]|2007l for the def¬ 
initions of these spectral states). There are two kinds of radio 
jets observed in BHXBs. One type is called continuous or com¬ 
pact jet, which is observed in the hard state and switched off in 
the high soft state. The continuous jets are compact and usually 
not easy to resolve with current radio f acilities, except in several 
cases, such as jets in GRS 1915+105 ( Dhawan et al. 1 12000 ) and 
Cyg X-l dStirling et al. Il200ll) . The main radiation component of 
the c ompact jet is thought to be self-abso rbed synchrotron emis¬ 
sion dMarkoff et al. l200ll : iYuan et al.l2005l) . and usually shows flat 
or slightly inverted spectrum in radio band. A correlation between 
the radio flux and the X-ray flux L^ad cx Ly 0,7 ) has been found 
in black hole X-ray binaries dHannikainen et all 


Corbel et al .1 


1200(1120031) and is thought to be universal 1 Gallo et al.l2003 ). How¬ 
ever, some outliers and an additional track have been found in re- 
cent years with much different correlation index (e.g. ISoleri et al.l 
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1 201)1 [Coriat et al. ||20 1 ll: iRatti et al.ll2012h . 


The other type is called episodic or ballistic jet, which 
displays intense radio flares during the intermediate state dur¬ 
ing the transitio n betw een the hard state and the high soft 
state t lFender. Belloni & Gallo I 120041) . At this time, the con¬ 
tinuous jet has already faded. The radio flares are associated 
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with plasmoids ejections, and sometimes can be resolved by 
radio observations, such as the jets seen in G RS 1915+105 
< Mirabel & Rodriguez 1 11994 IFender et alj Il999h . GRO J1655- 
40 lHjelhnm&Ru]3en 19951 : Tingavet al.1 199 31. XTE J1550- 
564 dCorbel et al.ll2002h. H 1743-322 dMiller-Jones et al.ll2012h 
and XTE J1752-223 jYang et aljEoIoL 1201 lh . The radio emis¬ 
sion of such episodic jet originates from an optically thin syn¬ 
chrotron emission component and the Lorentz factor of the ejected 
plasmoids is usually much larger than that of the continuous 
jet, and someti mes these episodic jets appea r superluminal (e.g . 


GRS 1915+105. IMirabel & Rodriguez 11 1994 IFender et all 


GRO J1655-40, Hjellmin & Runen 19951 : Hingav et al J 


1999 


1995; 


XTE J1748-288. iBrocksopp et a" 2007 ). The emission from the 
episodic j et is probably more polarized than the continuous jet 
jFender & Belloni 12004)) . 

Several models have been propos ed for powering continuous 
and episodic jets, from black hole spin iBlandfor d & Znajekl 1977b 
to diverse processes in t he accretion flow l|Falcke_&J3iennann] 
1995; Snad aet al J [200 1: Markoff et all 1200 ll: lYuan et alJ l200l 


200£). Whether black hole spin is related to the jet power of both 


continuou s or ep isodic j et has b ecome a central topic. For contin¬ 
uous jets. [Fender, Gallo & Russel] J2010h found no single relation 
between the black hole spin and the jet power. This means that 
the compact co nti nuous jets a re not powered by black hole spin. 
INaravan & McClintockl 120121) pointed out that the compact con¬ 
tinuous jet in the hard s tate was la unche d fa r away from the ra¬ 
dius (~ 10-100GM/c 2 . lM~arkoff et al]|2005l ) where the effect of 
the black hole spin is not important, so the spin has no signif¬ 
icant effect on the continuous jet power. However, the episodic 
jet was thought to be launched within several gravitational radii 
and_spin may play an important role in powering the episodic 
iet. lNaravan & McClintockl l2012h studied the relation between the 
black hole spin estimated via X-ray continuum-fitting method and 
the episodic jet power, which was defined as the peak luminosity of 
the radio flare at 5 GHz. They claimed that there is a positive corre¬ 
lation between the black hole spin and the power of the episodic jet 
in five black hole X-ray binaries, providing strong evidence that the 
episodi c jets were p owered by the spin of th e black hole. Lat er on , 
ISteiner et alj d2013l) followed the work of INaravan & McClintockl 
1 20121) bv~adding one more black hole X-ray binar y in the data 
and further strengthened this conclusion. However, [Russell et al .1 
J2013I ") reached a different conclusion. They have estimated the jet 
power with two different ways, one way was to estimate the peak 
jet power by taking the peak radi o luminosity of the radio flare, the 
same as INaravan & McClintockl (I2012I) . the other way was to esti¬ 
mate the minimum total energy required to energize a synchrotron 
flare. They found there was no relation between the jet power esti¬ 
mated as above and the black hole spin estimated via either X-ray 
continuum method or X-ray reflection method. Their results imply 
that the black hole spin does no t play an important role in powering 
the episodic jet. However, iMcClin tock et alj d2013h pointed out the 
data selection in lRussell et al J1 201 3l) was not appropriate because 
they neglected beaming effect. 

For discrete/episodic ejections seen in earlier observations of 
black hole binaries mostly in the transition phase (but in some 
cases the radio peak may be associated with X-ray low/hard 
state, see iFenden <20011) ). a significant correlation between the 
peak radio flux and the p e ak soft X - ray lu minosity was found 
( IFender & Kuulkersl 1200 11 ). lYu et al. <2004 0 pointed out there 
should be a correlation between the peak hard X-ray flux (i.e., hard- 
to-soft transition flux in most cases) and the peak radio flux in the 
hard state because there is a positive correlation between the hard 


X-ray peak flux and soft X-ray peak flux for several low mass X- 
ray binary transients. Such a correlation is expected if there is a 
universal correla tion between the radio flux and the hard X-ray flux 
in the hard state <Gallo et alj|2003l) . More studies on single sources 
has shown that there is a correlation between the transition lumi¬ 
nosity and the peak luminosity of the following soft state in single 
sources such as Aquila X—1 and GX 339—4 <Yu et al J1 2007a]fibh . 
Following that, a systematic study of bright persistent and tran¬ 
sient black hole and neutron star X-ray binari es with data from the 
RXTEI ASM and the Swift/BAT b^f Yu & Yard d2009l) demonstrated 
that the correlation found bv l Yu et al J <2004b holds for both flares in 
persistent sources and outbursts in transient sources. Furthermore, 
a correlation between the rate-of-change of the luminosity around 
the hard-to-soft state transition and the luminosity of the hard-to- 
soft state transition was found as well. The series of studies indicate 
that non-stationary accretion, represented by rate-of-change of the 
mass accretion rate instead of the mass accretion rate alone, plays a 
dominant role in the spectral state transitions and in producing the 
outburst flux amplitudes. Motivated by the series of works, it is very 
important to study whether non-stationary accretion, characterized 
by large rate-of-change in the mass accretion rate, plays important 
role in generating episodic jet in black hole X-ray binaries during 
the transition from the hard state to the soft state. 


2 OBSERVATIONS AND DATA ANALYSIS 

We intended to obtain the peak jet power and the peak luminos¬ 
ity of the soft state in the rising phase of any black hole transient 
outburst that has been observed by the RXTEI ASM. We adopted 
the peak luminosity of the r adio jet at 5 GHz as the peak p ower 
of the episodic jet following INaravan & McClintockl <2013) . The 
jet power is estimated by using the following formula: -Pjet — 
D 2 {vS„) max ,5GHz /M. This definition is actually corresponding 
to the radio jet luminosity L ra dio which is proportional to the jet 
power Qj et ( £, is in the range 1.22-1.5 8), when Qj e t is a co nstant 
fraction of accretion power Q aC c (see iFalcke & Biermamil 1 19951 : 


iHeinz & Sunvaevl [2003 Coriat£tal ii 


20UJ). We selected most 
of the sour ces from INara van & McClintock d2012l) . ISteiner et al] 
<20131) and iRussell et al.l <2013l) . except GRS 1915+105 and Cyg 


X-l. The X-ray observations of GRS 1915+105 allows us to search 
for the soft state which followed radio flares. But we found that 
no soft states can be associated with the radio flares in 1997, 2006 
and 2010, based on the hardnes s ratio between RXTEI ASM and 
Swift/BXT (see the method in I Yu & Van[|2009l) or the hardness ra¬ 
tio calculated with the RXTEI ASM 3-band data. For the 1994 radio 
flare we did not have good soft X-ray coverage, so do not allow us 
to determine the time range of the soft state. So we excluded GRS 
1915+105 in our plot, but we will discuss this source in §4. We also 
exclude Cyg X-l, because the sourc e didn’t enter the soft state after 
the radio flare on 2004 Feb 20l lFender et al.ll2006l) . 

In addition, we have also analyzed radio observations of the 
black hole transient GX 339-4 and found evidence of an episodic 
jet in the black hole transient GX 339-4 during the intermediate 
state of the 2010-2011 outburst, which has never been reported, al¬ 
though the measurements of the ra dio fluxes during the hard state of 
the outburst has been studied by lCorbel et al .1 <2013l) . We processed 
the radio observations of GX 339-4 observed by Australia Tele¬ 
scope Compact Array (ATCA) during its 2010-2011 outburst in 
order to obtain the measurement of the peak flux of the episodic jet. 
The ATCA synthesis telescope consists of six 22-m antennas along 
the east-west direction. It was upgraded with the new Compact Ar¬ 
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ray Broadband Backend (CABB) system dWilson et alj|201 ll) . and 
the maximum bandwidth had been increased from the 128 MHz 
(in each of two IF bands) to 2 GHz. This was coupled with higher 
level data sampling, and has improved the continuum sensitivity 
significantly. The observations of GX 339-4 were carried out at the 
frequencies of 5.5 and 9.0 GHz simultaneously. PKS 1934-638 was 
used as the flux and bandpass calibrator, and the antenna’s gain and 
phase calibration, as well as the polarization leakage, were carried 
out by nearby PKS 1740 -517. We used MIRIAD software package 
dSault & Killeen I H998I) to analyze the ATCA data. Our analysis 
of the observations corresponding to the hard state has reproduced 
the published results by Gorbel et alj d2013i) . We then focused on 
the radio data collected during its hard-to-soft transition in order to 
study the episodic jet. We found an optically-thin radio flare dur¬ 
ing the hard-to-soft state transition. The radio fluxes were averaged 
every 5 minutes for both bandpasses. We found that the radio flux 
increased from ~ 17 mjy to its peak flux at ~ 23.5 mjy at 5.5 GHz 
in less than 4 hours. The radio emission was optical thin during 
the entire observation. Figure 1 shows the ACTA light curve of the 
radio flare event on April 27, 2010, and the averaged flux on that 
day is shown on top right panel of Fig. 2 by the symbol of ma¬ 
genta star. The peak radio flux of the episodic jet in the outburst 
of 2010 is less than half of the pe ak r adio flux of a similar event 
in the 2002 outburst iGallo et a[]|2004h . The two outbursts clearly 
show the proportionality that brighter outbursts are associated with 
larger peak power of the episodic jets, implying the peak power 
of the episodic jet is primarily determined, if not entirely, by the 
accretion power that generates the X-ray outbursts. 

In order to compare the peak power of episodic jet and the 
peak luminosity of the soft state, we measured the corresponding 
peak fluxes of the soft states in the long-term light curves obtained 
with the RXTEI ASM (2-12 keV) for all the outbursts in black hole 
X-ray binaries with the peak episodic jet power of black hole tran¬ 
sients reported before. Figure 2 shows the daily averaged ASM 
light curves of the outbursts in nine black hole X-ray binaries in 
which both measurements of the episodic jet power and the soft 
X-ray peak are available. Data corresponding to the soft state are 
marked in red based on previous reports. In half of these reports, 
the d efinition of the soft state followed iRemilliad & McClintockl 
(2006 ). but in other reports earlier on the sources GX 339-4, H 
1743-322, XTE 11748-288, XTE J1720-318 and XTE J1652-453 
the definition of the soft state co uld be inconsistent with the defini¬ 
tion of lRemilliad & McClintockl J2006h . For the purpose to achieve 
consistent results, we have specifically studied the power density 
spectra obtained with the RXTEI PCA observations of those sources 
and measured the total root-mean-square (rms) power integrated 
over 0.1 to 10 Hz, in order to follow the same definition. The time 
ranges of the soft states we derived from those RXTE observations 
are almost the same as those reported w e applied the rms thresh- 
old fo r the soft state as less than 0.075 ( IRemilliad & McClintockl 
120061) and the differences are all within several days. We found 
these will not bring significant difference in our measurements of 
the peak flux of the soft state we selected based on the soft state 
regime identified in previous reports, except for XTE J1720-318, 
of which the time for the source to reach the peak flux of the soft 
state was two days later than the one from the previous report, but 
the flux difference can still be ignored due to larger uncertainties in 
the mass and distance. So the peak fluxes we measured are good es¬ 
timates of the peak fluxes of the soft states. The solid arrows mark 
the flux peaks of the soft states. There was only one exception - 
we excluded the isolated ASM flux peak in the soft state for XTE 


11752-223 in orde r to avoid f alse measurement of the peak flux, as 
discussed in dYu & Yanll2009l) . 

In order to accurately estimate the X-ray luminosity , we used 
the RXTEI ASM data in the energy bands 1.3~3.0 keV, 3.0~5.0 
keV and 5.0^12.2 keV. The Crab spectrum in 0.2-2 keV band, 2- 
10 keV band and 10-50 keV band can be described by a power law 
with photon indices of -2.02, -2.07 and -2.12, and no rmalizations of 
8.95, 8.26 and 9.42, respectively dKirsch et al.l200d ). We converted 
the count rates of the three bands of RXTEI ASM (in Crab unit) into 
luminosities in 2-12 keV by approximating these black hole X-ray 
binaries have simi lar spectral shape an d hydrogen absorptions to 
those of the Crab, lin’t Zand et al.1 ( 120071 ) estimated that, if the pho¬ 
ton index of the actual source differs to the Crab index by 1 or 
Nh is up to 10 times larger, the 2-12 keV energy flux differs by at 
most 30%. We have also simulated a disk-blackbody spectral com¬ 
ponent at 2 keV to see if we can use a Crab power-law spectrum 
to fit the data and to constrain the X-ray flux approximately in the 
three individual ASM bands with XSPEC. We found that the devi¬ 
ations are 30%, 4% and 2% in the 1.3-3.0, 3.0-5.0 and 5.0-12.2 
keV bands, respectively, while most of the X-ray flux comes from 
the upper two energy bands. So our estimates of the peak fluxes 
of the soft states are accurate (and consistent) enough for the cur¬ 
rent study. The corresponding peak luminosities (in Eddington unit, 
L^dd = 1.3 x 10‘ iS (M/ Mg)) of the high soft states of the black 
hole X-ray binaries are listed in Table 1. 

We have also managed to study the relation between the peak 
power of the radio episodic jet and the rate-of-increase of the X- 
ray luminosity during the rising phase in the hard state just before 
the spectral transition from the hard state to the intermediate state. 
We made use of the RXTEI ASM 3-band data to calculate the X-ray 
flux over the time window from the 30% of maximum ASM flux to 
the ASM maximum flux in the hard state and fit the curves with a 
linear relation, L = b + k x [t — 50000). Here, L is the luminosity 
in 2-12 keV, t represents time in unit of Modified lulian Day, and k 
represents the rate-of-change of the luminosity. The measurements 
of the rate-of-change of the X-ray luminosity are also shown in 
Table 1. In Fig. 2, the data used in the fits are shown as green circle, 
and the over-plotted blue solid lines correspond to the best-fits of 
the slopes, but here the slopes represent the rate-of-change of the 
X-ray flux. 


3 RESULTS 

3.1 Correlation between the episodic jet power and the peak 
luminosity of the soft state 

3.1.1 Galactic black hole transients 

Figure 3 shows the relation between the power of the episodic jet 
and the peak luminosity (in unit of Eddington luminosity) of the 
high soft state for all the sources in which an episodic jet power 
and the corresponding peak luminosity of the soft state can be es¬ 
timated with observations. Among the sources, we have relatively 
good estimates of the distances and black hole masses for H 1743- 
322, XTE J1550-564, XTE 11752-223 and GX 339-4 (The corre¬ 
sponding references are shown in Table 1.), which are plotted in 
black. Note that for XTE 11859+226, we only have the lower limit 
of black hole mass, we also plotted in black since the uncertainty in 
the mass would be within a factor of two. The distances and masses 
are not well known for XTE J1652-453, XTE J1748-288 and XTE 
J1720-318. The corresponding data are therefore plotted in green to 
distinguish them from good measurements. In the plot, we have set 
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the distance to 8 kpc and the mass to 10 solar masses for both XTE 
J1652-453 and XTE J1748-288 in the plot. For XTE J1720-318, 
our knowledge of source distance and mass is even worse, so we 
set the mass as 10 solar masses and the distance to be in the range 
from 3 kpc to 10 kpc, and the corresponding data are linked with a 
green dashed line in Fig. 3. Most of the es timates o f t he mass a nd 
distances were taken from previous wor ks l Naravan & McClintockl 


120121 : ISteiner et alj|2013l : iRussell et ai .1120 1 3I) . In order to increase 

our sample, we also looked into outbursts reported in A 0620-00, 
GS 1124-68, GS 2000+25 and GRO J1655-40 before the RXTE 
era. W e collected the peak fluxes in 0.4-10 keV from IChen et al.1 
d 1997h (note that whether the source are in the soft state is not 
clear). The corresponding luminosities in Eddington unit were then 
estimated with some updated source distances and masses (see Ta¬ 
ble 1). For the outburst of GRO J1655-40 in 1994, the peak flux 
in 0.4-10 keV was extr apolated from the reported peak flux in 20- 
100 keV dchen et al.ll997h . the source should have been in the hard 
state when the hard X-ray flux reached its peak. So the X-ray lumi¬ 
nosity in the 0.4—10 keV we estimated should be taken as the lower 
limit of the peak flux of the following soft state. The radio activi¬ 
ties of A 0620-00, GS 1124-68 and GS 2000+25 before the RXTE 
era were not well covered, which does not allow us to identify the 
radio flux peaks, therefore the measurements of the jet power also 
correspond to the lower limits. These samples are also plotted as 
green symbols in Fig. 3 and Fig. 4. 

When we only considered the sources with relatively good es¬ 
timates of distances and masses shown in black, which includes 
GX 339-4, H1743-322, XTE J1550-564 and XTE J1752-223, and 
XTE J1859+226, we found a strong correlation between the peak 
power of the episodic jet and the corresponding peak luminosity 
of the soft state. The Spearman rank correlation coefficient was 
0.86, with the chance possibility of 0.01. We also calculated the 
mean Spearman rank correlation coefficient and its chance pos¬ 
sibility when we considered the uncertainties of the peak X-ray 
luminosity of soft state and jet power using Bootstrap method. 
We found that the mean Spearman rank correlation coefficient is 
0.87±0.09, and the possibility is 0.02±0.03, implying a strong cor¬ 
relation between the peak power of the episodic jet and the peak 
luminosity of the soft state. We fit the data established for the 
seven outbursts of the above five sources with the linear function 
log Pjet = A + B x log L so ft,peak, and obtained the parameters 
A = 2.2±0.3 and B = 1.6±0.2. The best-fit model is over-plotted 
as a straight line in Figure 3. 

We should point out that the results shown in Figure 3 include 
the uncertainties estimated due to unknown distances or masses of 
black hole X-ray binaries. Notice that the green data points are sys¬ 
tematically below those data shown in black except XTE J1748-288 
and GRO J1655-40. This is consistent with the idea that the cover¬ 
age of the outbursts with radio observations which correspond to 
those green data points was much worse than the coverage with the 
X-ray observations due to the availability of X-ray monitors such 
as the RXTE/ ASM. Most of the radio flux peak corresponding to 
the episodic jet in those sources marked in green were probably 
missed due to sparse radio coverage. If we calculate the Spearman 
rank correlation for all the black hole X-ray binaries except XTE 
J1652-453 (which is put only a lower limit of the jet power) and 
set the distance of XTE J1720-318 to 6.5 kpc (taking as the me¬ 
dian of 3 kpc and 10 kpc - the range of previous measurements), 
we found the correlation coefficient is 0.47±0.07 and the chance 
possibility is 0.18±0.07. We found the source 4U 1543-475 has a 
significant impact on the result since the coverage of radio obser¬ 
vations of this source was bad and the peak radio flux was probably 


missed (we will discuss in next section). If we exclude this source, 
the spearman correlation coefficient would be 0.70±0.06 and the 
chance possibility would be 0.04±0.03. 

3.1.2 Extension to ultra-luminous X-ray sources in nearby 
galaxies 

The correlation between the power of the episodic jet and the corre¬ 
sponding peak luminosity of the soft state (or the rate-of-change of 
the X-ray luminosity) spans more than an order of magnitude in the 
parameter space without showing a luminosity saturation nor a cut¬ 
off at either ends of the X-ray luminosity range. This characteristic 
suggests that more luminous or dimmer episodic jet emission can 
be produced during more brighter or dimmer outbursts respectively, 
which hints that episodic jet production is somehow related to the 
occurrence of outbursts or specifically the rate-of-change in the 
mass accretion rate. This is similar to the conclusion obtained from 
the empirical relation between the hard-to-soft transition lumi nos¬ 
ity and the peak luminosity of the following soft state dYu & Yanl 
120091 ). Both empirical correlations probably extend to much higher 
luminosities such as those of the ultra-luminous X-ray sources 
(ULXs) in nearby galaxies l lYu & Yanl2009l) . On the other hand, as 
for the outbursts in the low mass X-ray binary (LMXB) transients, 
no saturation of the peak luminosity, up to the Eddington luminos¬ 
ity, is seen in the relation between the X-ray outbur st pe ak luminos¬ 
ity and the orbital period in the Galactic transients l lWu et al.l2010l) . 
These empirical relations suggested that some stellar mass compact 
objects can turn into ultra-luminous X-ray sourc es during brighter 
outbursts of similar duration or s horter outbursts 1 Yu & Yanll2009l) . 
which has been indeed observed dMiddleton et al.ll2013f ). 

The correlation between the peak power of episodic jet and the 
corresponding peak luminosity of the soft state in Galactic sources 
suggests that the investigation of the relation between the episodic 
jet power and the peak luminosity of the soft state in transient ULXs 
would tell us about the nature of the ULXs as compared with the 
Galactic black hole transients; whether ULXs contains intermedi¬ 
ate mass black holes or just extreme versions of the Galactic mi¬ 
croquasars. In order to address this question, we have investigated 
two transient ULXs in which radio emission from episodic jet has 
been detected during their outbursts. One of the ULXs is XMMU 
J004243.6+4125 Iff which wa s first d etected by XMM-Newton on 
15 January 2012 llHenze et al. 1120121} as an ultra-luminou s stellar- 
mass microquasar (~ 10 Mq, see Middleton et al.ll2013l l in M31, 
at a distance of only 0.78 Mpc away from us. Bright radio emis¬ 
sion was observed from this s ource w it h sim ilar behaviour to those 
Galactic black hole binaries l lMiddleton et alj|2013l) . The peak ra¬ 
dio flux of the first VLA observation ever measured was around 
0.5 mJy at 5.26 GHz. During several additional radio flares for this 
source the peak radio flux at 15 GHz was found around 1 mJy. 
So the peak radio flux at 5 GHz can be estimated roughly in the 
range 0.5-1 mJy. The peak X-ray luminosity of the soft state for 
the M31 ULX in 0.3-10 keV was around 1.26 xlO 39 erg s _1 
dMiddleton et al.ll2013h . We plotted this source as magenta open 
squares in Fig. 4. It is obvious that this source falls on the relation 
between the peak power of the episodic jet and the peak luminosity 
of the soft state established in the Galactic black hole transients as 
well. This implies that the correlation seems to hold up to around 
the Eddington luminosities for stellar mass black hole X-ray bina¬ 
ries. 

The other ULX detected in radio is HLX-1. HLX-1 was first 
observed by XMM-Newton on 23 November 2004 in the edge- 
on spiral galaxy ESO 243-49 about 95 Mpc away from us. The 
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measured maximum X-ray luminosity was about 1.1 x 10 42 erg 
s' 1 , which exceeds the Eddington luminosity of a stellar mass 
black hole by about three orders of magnitude. Radio observations 
have been performed with the Australia Telescope Compact Array 
(ATCA) during two spectral state tra nsitio n s of HLX-1, and radio 
flares were detected in this source ( jWebb et al.l2012h . The peak flux 
density of the radio flares for the combined 5 and 9 GHz data was 
63 ± 18 /rjy. This value should be the lower limit of the flux den¬ 
sity at 5 GHz given that similar radio flares observed in Galactic 
black hole systems are usually optical thin. We therefore set this 
value as the flux density at 5 GHz which should not deviate from 
the actual value too much. In order to obtain the peak luminosity 
of the soft state for HLX-1, we took the mass estimate of a black 
hole as between ~ 9 x 10 3 Mq and ~ 9 x 10 4 Mq iWebb et~ahl 
1 201 j) . We plotted these data of HLX-1 in Fig. 4 as blue open cir¬ 
cles with possible black hole mass of 9 x 10 4 Mq, 9 x 10 3 Mq and 
10 Mq, respectively (from left to right in sequence). Interestingly, 
these data fall on the correlation track formed from Galactic black 
hole binaries with good mass and distance measurements. The data 
point corresponding to the mass of 9 x 10 3 Mq is the closest point 
to the correlation track. This hints that HLX-1 is probably an in¬ 
termediate mass black hole if there exists a universal correlation 
between the peak power of the episodic jet and the peak luminos¬ 
ity of the soft state independent of black hole mass. However, we 
can not exclude that HLX-1 may contain a stellar mass black hole 
considering the 95% confidence interval of the fit in Fig. 4. 

3.2 Correlation between the peak power of episodic jet and 
the rate-of-change of the X-ray luminosity 

Figure 5 shows the relation between the power of the episodic jet 
and the rate-of-change of the X-ray luminosity around the hard- 
to-soft state transition measured from the ASM multi-band light 
curves. When we only considered the data with good estimates of 
X-ray flux shown in black, we found a very strong correlation. The 
Spearman rank correlation coefficient is 0.91±0.06, and the signif¬ 
icance is 0.008±0.013. If we take the distance and black hole mass 
of XTE J1748-288 as 8 kpc and 10 solar mass, respectively, the 
Spearman rank correlation coefficient is 0.98 and the chance for 
pure fluctuation is 3.3xl0~ 5 . This indicates there exists a strong 
correlation between the peak power of the episodic jet and the 
rate-of-change of the X-ray luminosity in black hole transient out¬ 
bursts. We fitted the data in black and that of XTE J1748-288 with 
a linear function of the form logPjet = A + B x \ogdL x /dt, 
here the unit of dL x /dt is Eddington luminosity per day, we got 
A = 2.4 ±0.2 and B = 0.9 ±0.1. If we excluded XTE J1748-288 
and only consider the data in black, we obtained A — 2.0 ± 0.4, 
and B = 0.7 ± 0.2. The best-fit model is over-plotted in Fig. 5 in 
red. 


4 DISCUSSION 

4.1 The accretion power and the episodic jet power 

Due to shorter accretion time scales and smaller distances, Galac¬ 
tic black hole binaries are the best targets to study episodic jets as 
compared with AGNs. Observations on time scales from days to 
years allow us to discriminate between persistent compact jets and 
episodic jets, which are usually seen in black hole hard state and in¬ 
termediate state, respectively. An important clue should be noticed 
in Fig. 3. In the black hole binaries H 1743—322 and GX 339—4, 


the episodic jet power and the peak luminosity of the soft state 
of the corresponding outburst have been measured for two distinct 
outbursts. In either sources, a larger power of the episodic jet dur¬ 
ing an outburst coincides with a larger peak luminosity of the soft 
state, which clearly indicates that the correlation is established in 
the accretion process among different outbursts when the black hole 
mass and spin should stays almost the same. Even if the total peak 
episodic jet power released might be somehow related to black hole 
spin, the majority of the episodic jet power should still come from 
the accretion power. In the current data which does not allow sub¬ 
sampling on black hole spin, the correlation established in outbursts 
of the same sources seems to be consistent with the overall correla¬ 
tion established in outbursts of our sample. This favors that a single 
physical process (parameter) drives the correlation. 

Since there is a correlation between the hard X-ray peak lumi¬ 
nosity (i.e., the luminosity of the hard-to-soft transition in outbursts 
with state transitions) and th e peak luminosity of the following soft 
state jYu et al .1120041 1 2007 a| [H; l Yu & Yanll20(M) . the correlation be¬ 
tween the peak episodic jet power and the peak luminosity of the 
soft state in this study implies a correlation between the peak lumi¬ 
nosity of the hard state (or the luminosity of the hard-to-soft state 
transition) and the peak power of the episodic jet. This seems to 
be supported by the apparent correlation between the power of the 
episodic jet and the accretion power measure d before episodic jet 
ejection seen in GRS 1915+105. IPunslv et alJl^Onl t found that the 
episodic jet power was strongly correlated with the intrinsic 1.2-50 
keV X-ray luminosity 0-4 hours before the episodic jet ejection, as 
well as the increase of the intrinsic X-ray luminosity at 1.2-50 keV 
during the hours before the ejection in GRS 1915+105. 

4.2 Episodic jets under non-stationary accretion 

Episodic jets are usually seen during spectral transitions between 
the hard state and the soft state when the source isjn the intermedi¬ 
ate state or steep power-law state dFender, Belloni & Gallo ll2004l ). 
It is known that there is a so-called ‘hyste resis’ effect of spectral 
state transitions during transient outburs ts dMivamoto et alj 1 19951 : 
lNowakll995 : Macca rone & Coppi 120031) . The dramatic hysteresis 
originate from the fact that the hard-to-soft spectral state transition 
occurs in a large luminosity range by u p to two order s of ma g- 
nitude in bright Galactic X-ray binaries dYu et al.l 1 20041 , 12007 al ibi: 
I Yu & Yaiill2009h . Systematic studies on the hard-to-soft transitions 
in both persistent and transient bright X-ray binaries show quan¬ 
titative evidence that non-stationary accretion, indicated by large 
rate-of-increase of the mass accretion rate during the rising phase 
of transient outbursts, plays an important role in determination of 
the luminosity at which the hard-to-soft transition occurs as well 
as the outburs t peak luminosity (or peak luminosity of the follow¬ 
ing soft statel dYu & Yanll2009r) . In this work, we have found that 
there is a correlation between the peak power of the episodic jet 
and the corresponding peak luminosity of the soft state, as well as a 
correlation between the peak power of the episodic jet and the cor¬ 
responding rate-of-increase of the X-ray luminosity in the 2-12 keV 
band during the rising phase of a transient outburst, indicating that 
non-stationary accretion, represented by the observed large rate- 
of-increase of the luminosity, appears to play an important role in 
generating the episodic jet and the determination of its peak power. 
We show that in multiple outbursts of a single source, i.e GX339-4 
and HI743-322, the larger the rate-of-increase of the X-ray lumi¬ 
nosity in the hard state is observed, the higher the peak power of 
the episodic jet will be (see Fig. 5). 

It is worth noting that the peak luminosity of the soft state and 
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the peak power of the episodic jet are not related to each other lin¬ 
early, thus the correlation cannot simply be explained as the likely 
connection between the instantaneous mass accretion rate and the 
mass loss rate into the jet. The relation between the peak luminos¬ 
ity of the soft state and the peak power of the episodic jet follows 
the relation Pj et ~ ±s 0 ft, P eak- Assuming the broadband spectrum 
of the episodic jet emission is of a power-law with an index of - 
0.5 all the way up to the optical band, the total jet luminosity is 
about 500 times the luminosity at 5 GHz, while the total X-ray lu¬ 
minosity in 0.1-100 keV is about 12 times the X-ray luminosity 
in 2-12 keV if we assume the photon index is 1.7. We found the 
prompt episodic jet luminosity will surpass the X-ray luminosity 
when the X-ray luminosity in 0.1 -100 keV is larger than nine times 
Eddington luminosity. This gives an order-of-magnitude estimate 
of the luminosity of the non-stationary regimes corresponding to 
such a possible episodic jet-dominated state. The luminosity limit 
also serves as a threshold to discriminate whether the episodic jet 
power would primarily come from an energy reservoir, e.g., the en¬ 
ergy from the spinning black hole by the Blandford-Znajek (BZ) 
effect or magnetic energy in the accretion flow, or instead from 
the instantaneous accretion power. Both the episodic jet-dominated 
state and the threshold are interesting predictions by the empirical 
relations in the parameter space not accessible now. 

There are several ways to understand the relation between the 
peak luminosity of the soft state and the peak power of the episodic 
jet, which follows the relation Pj e t ~ -i'aoft, P eak> when we rec¬ 
ognize the role of the rate-of-change of the mass accretion rate in 
interpreting the observations unaccounted for by the mass accretion 
rate alone. One possibility is that, if the radiation efficiency of the 
mass loaded in the episodic jet does not change, then the ratio of the 
mass loss rate in the jet to the mass accretion rate in the accretion 
flow should increase with outburst amplitude and correlate with the 
rate-of-change of the mass accretion rate. The typical fraction of 
the mass lose rate into jet in those outbursts we studied as compared 
with the mass accreti on rate is probably only a few percent, such as 
0.5% (see Yuan et al . 200gl. Then this mechanism can at most con¬ 
tribute to a parameter space spanning by ~ 200 times, from 0.5% 
to 100% of the mass accretion rate. The other possibility is that the 
ratio of the mass loss rate in the jet to the mass accretion rate in 
the accretion flow remains roughly constant but the radiation effi¬ 
ciency of the episodic jet increases with peak luminosity of the soft 
state and correlates with the rate-of-change of the mass accretion 
rate. The limiting factor of this mechanism would be on the order 
of 100 times as well, since the radiative efficienc y of jet is proba¬ 
bly on the order of a few p ercent, such as 0.05 l lFender et all 200 ll : 
lYuan et ail 20051) or 0.003 l lMalzac et all 20041) . Yet another possi¬ 
bility is that the enhanced episodic jet power actually comes from 
the black hole through B-Z effect, which requires the extraction of 
black hole energy more efficient for sources with brighter outbursts 
and for sources with higher spin. If this is true, we would see sys¬ 
tematic shifts of peak power of episodic jet among black holes of 
different spins, forming parallel tracks of the relation between the 
peak power of episodic jet and peak luminosity of the soft state in 
Fig. 3. This is yet not seen in the current data, but an increase of 
the sample to several tens outbursts of sources with diverse spins 
would allow meaningful constraints. 

Spinning black hole as compared with non-spinning black 
holes would have higher radiation efficiency for the emission of the 
the episodic jet in the radio band as well as for the accretion flow 
in X-rays by similar factors. This suggests that if we plot multiple 
outbursts of individual black hole binaries in Fig. 3 or Fig. 4, black 
hole binaries with different black hole spin would not bring sig¬ 


nificant scattering around the overall empirical relation we found, 
but only bring shifts up and down along the empirical relation. If 
the power of the episodic jet primarily comes from the energy ex¬ 
traction from the spinning black hole, we would observe diverse 
correlations between the peak episodic jet power and the peak lu¬ 
minosity of the soft state among sources with significantly different 
black hole spins. This in turn suggests that towards the higher lu¬ 
minosity end in the relation shown in Figure 3, a larger scatter of 
observations data should exist if the episodic jet power is strongly 
dependent on the black hole spin. Current observations do not show 
such a sign. 


4.3 Episodic jets: early prediction of its peak radio flux 

lYu & Ya5 J20091 ) found the relation between the luminosity of 
hard-to-soft state transition and the peak luminosity of the follow¬ 
ing soft state in Eddington unit is in the form log L so ft,peak = 
(0.4±0.1) ± (0.9±0.1) log Ltr, where Ltr represents the hard-to- 
soft state transition luminosity in 15-50 keV in unit of Eddington 
luminosity. So the relation between the peak power of the episodic 
jet and the luminosity of hard-to-soft state transition is derived as 
log Pjet = (2.9± 0.4) ± (1.5 ± 0.3) log Ltr-These empirical rela¬ 
tions derived above are useful. We are able to predict the peak flux 
of the episodic jet based on the luminosity of the hard-to-soft state 
transition and the rate-of-increase of the X-ray luminosity before 
an episodic jet is launched. For example, when the luminosity of 
the hard-to-soft state transition is only 0.001 L^dd and the black 
hole of 10 solar masses is at a distance of 8 kpc away from us, 
we can roughly estimate that the peak flux of the episodic jet at 5 
GHz is ~0.6 mJy. If the hard-to-soft transition occurs at 0.1 ±Edd, 
then the peak flux of the episodic jet at 5 GHz would be ~ 0.7 Jy 
(see lYu & Yanll2009l) . In some extreme cases when the hard-to-soft 
transition occurs at nearly L/Edd and a distance of 0.78 Mpc (M31) 
or 3.25 Mpc (M82) for ultraluminous X-ray sources, then the peak 
flux of the episodic jet at 5 GHz would be ~2.3 mJy or ~0.1 mJy, 
respectively. Normally episodic jet flares in microquasars are on the 
time scale of hours, thus sensitive wide field-of-view radio facilities 
such as Murchison Widefield Array (MWA) and Square Kilometer 
Array (SKA) would be able to identify such extreme events in black 
hole X-ray binaries in nearby galaxies and efficiently probe the na¬ 
ture of ultra-luminous X-ray sources in nearby galaxies. 


5 SUMMARY 

In this paper, we have collected published radio data of episodic 
jets observed in nine black hole X-ray binaries which were covered 
simultaneously by the X-ray monitoring with the RXTE/ASM. In 
addition, we have included our own detection of an episodic jet in 
the ATCA radio observations of GX 339-4 in its 2010 outburst dur¬ 
ing the hard-to-soft state transition. Using these data, we found a 
strong correlation between the peak power of the episodic jet and 
the peak X-ray luminosity of the soft state among the outbursts. 
Data from multiple outbursts of individual sources and data among 
all the sources seem to agree with a single relation, indicating the 
power of the episodic jet is primarily accretion-driven. It is surpris¬ 
ing that two transient UFXs found in nearby galaxies in which radio 
flares were detected during spectral state transitions also fall on the 
same track. We also found a correlation between the peak episodic 
jet power and the rate-of-increase of the X-ray luminosity among 
the sources and among outbursts. These discoveries imply that the 
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power of the episodic jet is primarily driven by non-stationary ac¬ 
cretion during the rising phase of transient outbursts, the scale of 
which varies from outburst to outburst and from source to source. 

Our results not only confirm the previous discovery of a 
positive correlation between the peak X-ray flux of X-ray out- 
burst and the peak radio flux density for black hole binaries 
( [Fender & Kuulkersjl200lh . but also reveal the importance of rate- 
of-change in the mass accretion rate to the episodic jet power. In 
this paper we intend to focus on the origin of the power of the 
episodic jet in black hole transients. We found two correlations, 
namely the correlation between the peak power of the episodic jet 
and the peak luminosity of the soft state, and the correlation be¬ 
tween the peak power of the episodic jet and the rate-of-increase 
of the X-ray luminosity around the hard-to-soft state. Our results 
provide strong evidence that non-stationary accretion is the driving 
source of the episodic jets seen in black hole transients. Models 
about jet production should consider the influence of the apparent 
rate-of-change of the mass accretion rate shown in the change of 
the X-ray flux. 


ACKNOWLEDGMENTS 

We thank the anonymous referee’s comments, and we also thank 
Jean-Pirre Lasota, Feng Yuan and Qinwen Wu for their good 
suggestions and comments. We thank the RXTE and the Swift 
Guest Observer Facilities at NASA Goddard Space Flight Cen¬ 
ter for providing the RXTE /ASM products and the Swift/BXT 
transient monitoring results. This work was supported in part by 
the National Natural Science Foundation of China (245681001, 
11333005, 11073043, and 11350110498) and by Strategic Prior¬ 
ity Research Program ’’The Emergence of Cosmological Struc¬ 
tures” under Grant No. XDB09000000 and the XTP project un¬ 
der Grant No. XDA04060604, by the Shanghai Astronomical Ob¬ 
servatory Key Project. This work has made use of the data ob¬ 
tained through the High Energy Astrophysics Science Archive Re¬ 
search Center Online Service, provided by the NASA/Goddard 
Space Flight Center. This paper also includes analysis of archival 
data obtained through the Australia Telescope Online Archive 
(http://atoa.atnf.csiro.au I. 


REFERENCES 

Ball, L., Kesteven, M. J., Campbell-Wilson, D., Turtle, A. J., 
Hjellming, R. M. 1995, MNRAS, 273, 722 
Barret, D., McClintock, J. E., & Grindlay, J. E. 1996, ApJ, 473, 
963 

Blandford, R. D„ & Znajek, R. L. 1977, MNRAS, 179, 433 
Brenneman, L. W. et al. 2011, ApJ, 736, 103 
Brocksopp, C. et al. 2002, MNRAS, 331, 765 
Brocksopp, C. et al. 2005, MNRAS, 356, 125 
Brocksopp, C., Corbel, S., Tzioumis, A., Broderick, I. W., Ro¬ 
driguez, J., Yang, I., Fender, R. P., & Paragi, Z. 2013, MNRAS, 
432, 931 

Brocksopp, C., Miller-Jones, I. C. A., Fender, R. P., & Stappers, 
B. W. 2007, MNRAS, 378, 1111 

Calvelo, D. E., Tzioumis, T., Corbel, S., Brocksopp, C., Fender, 
R. P. 2009, Astron. Telegram, 2135 
Cantrell, A. G. et al. 2010, Apl, 710, 1127 
Chaty, S„ & Bessolaz. N. 2006, A&A, 455, 639 
Chen, W„ Shrader, C. R„ & Livio, M. 1997, Apl, 491, 312 


Corbel, S., Broderick, J., Brocksopp, C., Tzioumious, T., Fender, 
R. 2010, Astron. Telegram, 2525 

Corbel, S., Fender, R. P., Tzioumis, A. K., Nowak, M., McIntyre, 
V., Durouchoux, P, & Sood, R. 2000, A&A, 359, 251 
Corbel, S.. Fender, R. R, Tzioumis, A. K., Tomsick, J. A., Orosz, 
J. A., Miller, I. M., Wijnands, R., & Kaaret, P. 2002, Science, 
298, 196 

Corbel, S., Nowak, M. A., Fender, R. R, Tzioumis, A. K., & 
Markoff, S. 2003, A&A, 400, 1007 
Corbel, S., Coriat, M., Brocksopp, C., Tzioumis, A. K., Fender, 
R. P., Tomsick, J. A., Buxton, M. M., & Bailyn, C. D. 2013, 
MNRAS, 428, 2500 

Coriat, M., Corbel, S., Prat, L., Miller-Jones, J. C. A., Cseh, D., 
Tzioumis, A. K., Brocksopp, C., Rodriguez, J., Fender, R. P., & 
Sivakoff, G. R. 2011, MNRAS, 414, 677 
Corral-Santana, J. M., Casares, J., Shahbaz, T., Zurita, C., 
Martlnez-Pais, I. G., Rodrlguez-Gil, P. 2011, MNRAS, 413, L15 
Dhawan, V., Mirabel, I. F., & Rodrlgure, L. F. 2000, ApJ, 543, 373 
Dinner, T., Kalemci, E., Buxton, M. M., Bailyn, C. D., Tomsick, 
J. A„ Corbel, S. 2012,ApJ, 753, 55 
Done, C., Gierlinski, M., & Kubota, A. 2007, A&AR, 15, 1 
Falcke, H. & Biermann, L. 1995, A&A, 293, 665 
Fender, R. P. 2001, MNRAS, 322, 31 
Fender, R„ Belloni, T. 2004, ARA&A, 42, 317 
Fender, R„ Belloni, T„ Gallo, E. 2004, MNRAS, 355, 1105 
Fender, R„ Gallo, E„ Russell, D. 2010, MNRAS, 406, 1425 
Fender, R. P., Garrington, S. T., McKay, D. J., Muxlow, T. W. B., 
Pooley, G. G., Spencer, R. E., Stirling, A. M., & Waltman, E. B. 
1999, MNRAS, 304, 865 

Fender, R. R, & Kuulkers, E. 2001, MNRAS, 324, 923 
Fender, R. P. et al. 2001, MNRAS, 322, L23 
Fender, R. R, Stirling, A. M., Spencer, R. E., Brown, I., Pooley, 
G. G„ Muxlow, T. W. B„ Miller-Jones, J. C. A. 2006, MNRAS, 
369, 603 

Gallo, E„ Fender, R. P„ & Pooley, G. G. 2003, MNRAS, 344, 60 
Gallo, E., Corbel, S., Fender, R. P, Maccarone, T. J., & Tzioumis, 
A. K. 2004, MNRAS, 347, 52 

Gelino, D. M., Harrison, T. E., & McNamara, B. J., 2001, AJ, 122, 
971 

Greene,J., Bailyn, C. D., & Orosz, J. A. 2001, ApJ, 554, 1290 
Hannikainen, D. C., Hunstead, R. W., Campbell-Wilson, D., & 
Sood, R. K. 1998, A&A, 337, 460 
Hannikainen, D. C. et al. 2009, MNRAS, 397, 569 
Heinz, S. & Sunyaev, R. A. 2003, MNRAS, 343, L59 
Henze, M. et al. 2012, Astron. Telegram, 3890 
Hjellming. R. M., Calovini, T. A., Han, X. H., Cordova, F. A. 
1988, ApJ, 335, L75 

Hjellming. R. M., Rupen, M. P. 1995, Nature, 375, 464 
Homan, J., Buxton, M., Markoff, S., Bailyn, C., Nespoli, E., Bel¬ 
loni, T. 2005, ApJ, 624, 295 
Hynes, R. I. 2005, ApJ, 623, 1026 

in’t Zand, J. J. M., Jonker, P. G., Markwardt, C. B. 2007, A&A, 
465, 953 

King, A. L„ et al. 2013, ApJ, 771, 84 
Kirsch, M. G. et al. 2005, Proc.SPIE, 5898, 22 
Kording, E. G., Jester, S., & Fender, R. 2006, MNRAS, 372, 1366 
Kuulkers, E., Fender, R. R, Spencer, R. E., Davis, R. J. & Morison, 
I. 1999, MNRAS, 306, 919 

Liu, Q. Z., van Paradijs, J., van den Heuvel, E. P. J. 2007, A&A, 
469, 807 

Ma, R. Y. 2012, MNRAS, 423, L87 

Maccarone, T. J., & Coppi, P. S. 2003, MNRAS, 338, 189 


© 0000 RAS, MNRAS 000, 000-000 




8 Zhang & Yu 


Malzac, J., Merloni, A., & Fabian, A. 2004, MNRAS, 351, 253 
Markoff, S„ Falcke, H„ & Fender. R. 2001, A&A, 372, 25 
Markoff, S„ Nowak, M. A.; Wilms, J. 2005, ApJ, 635, 1203 
McClintock, J. E., Remillard, R. A., Rupen, M. P., Torres, M. A. 

P., Steeghs, D., Levine, A. M., Orosz, J. A. 2009, ApJ, 698, 1398 
McClintock, J. E., Narayan, R., Steiner, J. F. 2013, Space Science 
Reviews (arXiv: 1303.1583) 

Middleton, M. J. et al. 2013, Nature, 493, 187 
Miller-Jones, J. C. A. et al. 2012, MNRAS, 421, 468 
Mirabel, I. F., Rodriguez, L. F. 1994, Nature, 371, 46 
Miyamoto, S., Kitamoto, S., Hayashida, K., Egoshi, W. 1995, 
ApJ, 442, L13 

Moderski, R„ Sikora, M„ & Lasota, J. R 1998, MNRAS, 301, 142 
Motta, S., Belloni, T., Munoz Darias, T. 2010, Astron. Telegram, 
2545 (2010a) 

Motta, S., Belloni, T., Munoz-Darias, T.. Homan, J. 2010, Astron. 

Telegram, 2593 (2010b) 

Nakahira, S. et al. 2012, PASJ, 64, 13 

Narayan, R., Mahadevan, R., & Quataert, E. 1998, in Abramowicz 
M. A., Bjornsson G., Pringle J. E., eds, Theory of Black Hole Ac¬ 
cretion Disks, Cambridge Univ. Press, Cambridge, U. K. p.148 
Narayan, R„ & McClintock, J. 2012, MNRAS, 419, L69 
Nowak, M. A. 1995, PASP, 107, 1207 

Orosz, J. A., Polisensky, E. J., Bailyn, C. D., Tourtellotte, S. W., 
McClintock, J. E„ Remillard, R. A. 2002, BAAS, 34, 1124 
Orosz, J. A., Steiner, J. F., McClintock, J. E., Torres, M. A. P, 
Remillard, R. A., Bailyn, C. D„ Miller, J. M„ 2011, ApJ, 730, 
75 

Park, S. Q., et al. 2004, ApJ, 610, 378 
Punsly, B., & Rodriguez, J. 2013, ApJ, 764, 173 
Ratti, E. M. et al. 2012, MNRAS, 423, 2656 
Remillard, R. A. & McClintock, J. E. 2006, ARA&A, 44, 49 
Revnivtsev, M. G., Trudolyubov, S. P, Borozdin, K. N. 2000. MN¬ 
RAS, 312, 151 

Russell, D. M„ Gallo, E„ Fender, R. P. 2013, MNRAS, 431, 405 
Sault, R. J., & Killeen, N. E. B. 1998, The miriad User’s Guide. 

Sydney: Australia Telescope National Facility 
Shaposhnikov, N., Markwardt, C., Swank, J., & Krimm, H. 2010, 
ApJ, 723, 1817 

Shaposhnikov, N., & Titarchuk, L. 2009, ApJ, 699, 453 
Smith, D. M., Kalemci, E., Tomsick, J. A., Heindl, W. A. 2003, 
Astron. Telegram, 120 
Sobczak, G. J. et al. 2000, ApJ, 544, 993 
Soleri, P. et al. 2010, MNRAS, 406, 1471 
Spada, M., Ghisellini, G., Lazzati, D., & Celotti, A. 2001, MN¬ 
RAS, 325,1559 

Steiner, J. F., McClintock, J. E., & Narayan, R. 2013, ApJ, 762, 
104 

Steiner, J. F„ McClintock, J. E„ Reid, M. J., 2012, ApJ, 745, L7 
Stirling, A. M., Spencer, R. E., de la Force, C. J., Garrett, M. A., 
Fender, R. P„ & Ogley, R. N. 2001, MNRAS, 327, 1273 
Tang, J., Yu, W. F„ & Yan, Z. 2011, RAA, 11, 434 
Tingay, S. J. et al. 1995, Nature, 374, 141 
Webb, N. et al. 2012, Science, 337, 554 
Wilson, W. E et al. 2011, MNRAS, 416, 832 
Wilson, A. S., & Colbert, E. J. M. 1995, ApJ, 438, 62 
Wu, Y. X., Yu, W„ Li, T. P, Maccarone, T. J., Li, X. D. 2010, ApJ, 
718, 620 

Yamaoka et al., 2010, Astron. Telegram, 2380 
Yang, J., Brocksopp, C., Corbel, S., Paragi, Z., Tzioumis, T., & 
Fender, R. P 2010, MNRAS, 409, 64 


Yang, J., Paragi, Z., Corbel, S., Gurvits, L. I., Campbell, R. M., & 
Brocksopp, C. 2011, MNRAS, 418, 25 
Yu, W„ van der Klis, M„ Fender, R. 2004, ApJ, 611, L121 
Yu, W., Lamb, F. K., Fender, R. P, van der Klis, M. 2007, ApJ, 
663, 1309 

Yu, W„ Dolence, J. 2007, ApJ, 667, 1043 

Yu, W., & Yan, Z. 2009, ApJ, 701, 1940 

Yuan, F„ Cui, W., & Narayan, R. 2005, ApJ, 620, 905 

Yuan, F„ Lin, J., Wu, K, Ho, L. C. 2009, MNRAS, 395, 2183 

Zdziarski, A. A„ et al. 2004, MNRAS, 351, 791 

Zhang, H„ Yuan, F„ & Chaty, S. 2010, ApJ, 717, 929 


© 0000 RAS, MNRAS 000, 000-000 


Episodic jet power and peak luminosity of soft state 9 



012345678 
Time From Observation Start (Hour) 


Figure 1. The radio flare of GX 339-4 seen with the ATCA on 2010 April 27 (MJD 55 313) at 5.5 and 9.0 GHz (upper) and the corresponding spectral index 
a (lower, where S u oc 1 /“). The error bars are overlapped by the symbols. The radio fluxes are averaged over 5 minutes. 
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Figure 2. X-ray and radio observations of a sample of outbursts seen with the RXTEI ASM in the RXTE era. The vertical dashed line shows the time of the radio 
peak flux from the episodic jet. The blue stars and green filled squares represent flux measurements from radio observations , the magenta stars in to p middle 
panel represents the peak of the radio flare event of GX339-4 in 2010. The ra dio light curve of X TE J1859+226 can refer tolBrocksopp et a lj |2002|) . The re d 
filled circles represent data corresponding to the high soft stat e, for GX 339-4 jHomaneta lJ2005HSmith et al.l2003UM otta et all20Kp lhHpincer et all2012l> . 
HI743-322 jMcClintock et all2009tlMiller-Jones et all2012h. XTE J1752 -223 jNakahira et all2012h. XTE J1550-564 ISobczak et all200Ch . XTE J1859+226 
Remilliari & McClintocklkOOd) . XTE J1720-318 iBrocksopp et al.ll2005h . XTE J1748-288 jBrocksopp et all2007HRevnivtsev et all200Ch and 4U 1543-475 
Park et all20041) . Please note that the intermediate state could be brighter than the soft state in several black hole transient outbursts in 2-12 keV band, which 
cover both thermal disk component and power-law component. The data marked as green circles were fit to a straight line (shown in blue) to measure its slope 
and to derive the rate-of-increase of the X-ray luminosity in hard state. For 4U1543-475 and XTE J1652-453, no reliable data are available to measure the 
rate-of-change of X-ray luminosity. Since radio emission from the episodic jets comes from ejected blobs, in a few cases the radio flux peak is found when the 
source had entered the soft state. 
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Figure 3. The relation between the episodic jet power and the peak luminosity of the soft state in black hole transients. The black points correspond to sources 
have good estimates of mass, distances. X-ray and radio fluxes. The only exception is XTE J1859+226, which has only a lower limit of black hole mass. 
The green data points correspond to sources with a lot larger uncertainties in the estimates of their mass or distances. We show XTE J1720-318 with possible 
distances of 3, 6.5 and 10 kpc and the corresponding data points are connected with a green dashed line. Data corresponding to two outbursts of GX 339-4 
and H 1743-322 demonstrate the dominate role of accretion in producing radio jet power. For the sources A 0620-00, GS 1124-68, GS 2000+25 and GRO 
J1655-40 in which radio flares were reported before the RXTE era, the peak X-ray luminosity are calculated based on the extrapolations to the 0.4-10 keV 
band from hard state hard X-ray data, and the radio fluxes should correspond to the lower limits on the peak radio jet power due to very sparse coverage at 
early times. The red solid line are the best-fit for the black points in logarithm scale, and the red dashed lines are the 95% confidence intervals of the fit. The 
error bars of each data points are primarily from the uncertainties in the estimates of their distance and masses. 
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Figure 4. The relation between the episodic jet power and the peak luminosity of the soft state in black hole transients and a few ULXs. We add the data of 
HLX-1 and the ULX XMMU J004243.6+412519 in M31 to illustrate that the ULXs in nearby galaxies are also consistent with the relation found in Galactic 
black hole X-ray binaries. The blue circles which are connected with blue dotted line corresponds to different black hole mass estimates of HLX-1. From left 
to right, we took the black hole mass as 9 x 10 4 Mq, 9 x 10 3 Mq and 10 Mq. The magenta open squares correspond to radio measurements of 0.5 and 1 
mJy for the ULX in M31. The red solid and dashed lines are the same as that in Fig. 4. 
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Figure 5. The relation between the peak power of the episodic jet and the rate-of-change of the X-ray luminosity in the rising phase of the outbursts. The data 
shown in black correspond to sources with rather definite ranges of masses, distances, X-ray and radio fluxes, except for XTE J1859+226 which has only a 
lower limit of black hole mass. The green data points indicate sources a lot more larger uncertainty in the estimates of their masses and distances. In the plot, 
we set the distance to 6.5 kpc for XTE J1720-318. The red solid line are the linear fit for the data in black in logarithm scale, and the red dashed lines are the 
95% confidence intervals of the fit. 
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Table 1 . List of the sources and parameters used 


Source 

d 

(kpc) 

Mass 

(M 0 ) 

<S'5GHz 

(mJy) 

Lmax/LEdd 

(%) 

dL/dt 

(10 —2 L Bdd /Day) 

Year 

d 

References 

Mass S 5G Hz 

4U 1543-475 

7.5 3= 1.0 

9.4 3= 1.0 

17.5 

44.23=8.4 

— 

2002 

1 

1 

2 

GX 339-4 

83=2 

12.3±1.4 

~ 55 

8.63=2.7 

0.093=0.03 

2002 

3 

4 

5 




~24.5 

3.93=1.3 

0.063=0.02 

2010 




H 1743-322 

8.53=0.8 

13.33=3.2 

96.1 

8.43=2.3 

0.513=0.21 

2003 

6 

4 

7 




24.4 

2.63=0.7 

0.153=0.07 

2009 



8 

XTE J1550-564 

4.38±0.50 

9.103=0.61 

265.0 

8.53=1.3 

0.723=0.21 

1998 

9 

9 

10,11 

XTEJ1652-453(“) 

8 

10 

>0.95 

2.43=0.2 

- 

2009 

__a 

_a 

12 

XTEJ1720-318W 

3~10 

10 

4.7 

3.63=0. lM 

0.373=0.39 

2003 

13 

_b 

14 

XTEJ1748-288(“) 

8 

10 

500 

6.13=0.1 

3.693=0.11 

1998 

_a 

_a 

15 

XTE J1752-223 

3.5 3= 0.4 

9.6 3=0.9 

20 

1.43=0.2 

0.013=0.02 

2010 

16 

16 

17 

XTE J1859+226 

8±3 

>5.42 

100 

14.03=6.1 

2.353=1.78 

1999 

18 

19 

20 

GRO J1655-40 

3.23=0.5 

6.303=0.27 

2420^ 

>6.7 

- 

1994 

21 

22 

21 

A 0620-00 

1.06±0.12 

6.613=0.25 

>203< d > 

15.63=3.5 

- 

1975 

23 

23 

24 

GS 1124-68 

5.0 

6.03=1.5 

>17l(d) 

67.53=0.0 

- 

1991 

25 

26 

27 

GS 2000+25 

2.7±0.7 

7.23=1.7 

>4 g(d) 

18.13=9.4 

- 

1988 

28 

28 

29 


assume a distance of 8 kpc and a black hole mass of 10 Mq. assume a distance of 3 kpc and 10 kpc, respectively, and a black hole mass of 10 Mq. 

assume a distance of 6 .5 kpc. ^ measured in 0.4- 1 0 keV . 

(l lOrosz et a i] fc002l): (2lPark et al l <2004: (3lZdziarski et all fc004l): (4lShanoshniknv ft Titarchukl <20091) : (5lGallo et all <2004: ( 6 lSteiner et ail <20121) : 


l'7lMcClintock et alj ^20091): ( 8 )lMiller-Jones et alJ <20121): (9, 

fl2 lCalvelo et all J2009t): (13lChatv & Bessolaz N2006 ) ;(14. _ 

fl7lBrocksoDDet alJ ll2013h: d81HvnesTl2005h: (19 CorraPSanIana^H^"| 2 ^]]© 20 ^^o^^o ct al .| f 2002 |):l 21 jHjcllmin & Ruji 


OroszetaJ 2011): (lOlHannikainen et al.N2009t): (lllNaravan & McClintocklj2012l)~ 

jto^^TeuiiT 20(d); (15 (16 Shagoshnikovgtal j2O10j); 

N 19951): 


(22 lGreene et al.U200ll) :(23 lCantrell et al.l l2()lfl); (2 4lKuulkereetalTjl999l):(25jGelino et al J 120oT ): (26 lRemilliad &~McCUntoc3T20Mk (27 lBall et alj 

Jl995h : (28 lBarret et all il99ai :(29 lHiellming et al.l <198il . 
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